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A genuinely non-perturbative evaluation of higher-twist contributions to the structure functions of the nucleon, 
with all mixing effects and renormalon ambiguities taken care of, is presented. Higher-twist corrections turn out 
to be significant at moderate values of q^. 

1. INTRODUCTION 

The calculation of power corrections to the 
deep-inelastic structure functions of the nucleon 
is a long-standing problem. Recent phenomeno- 
logical studies [ |l| have shown that they account 
for more than 10% of the lower moments of F2 at 
q^ = 5 GeV . For the theoretical understanding 
of the structure functions, in particular their evo- 
lution as a function of q^ , it is thus important to 
gain quantitative control of these effects. 

The theoretical basis for the calculation is the 
operator product expansion (OPE): 
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where m distiguishes different operators with the 
same Lorentz symmetries. The Wilson coeffi- 
cients C are universal and depend only on q and 
the renormalization point /z, while all dependence 
on the target momenta is contained in the opera- 
tor matrix elements. The product of Wilson coef- 
ficients and matrix elements does not depend on 
/i, nor on the renormalization scheme. 



Besides the familiar operators of twist two, 
such as 

Cpi-'-Mn = '/'7{mi-C'm2 • • • ^M>.}'/' ~ traces, 
the OPE (|^) receives contributions from twist- 
four and higher operators, which give rise to 
power corrections, and which so far have been 
neglected in the theoretical analysis. A typical 
twist-four operator is 

Omi-m,x = ■07a'i-D/j2 ■ • --Da-Da • • ■-Dp„'0- 

In addition to the calculation of matrix ele- 
ments of higher-twist operators, the evaluation 
of power corrections requires a non-perturbative 
computation of the Wilson coefficients. This is 
necessary because of renormalon ambiguities and 
mixing effects, which tie Wilson coefficients and 
higher- twist matrix elements together [ 0] . 

In [ ^ we have developed a method that allows 
a lattice calculation of the Wilson coefficients. 
With this method we can compute the r.h.s. of 
P) entirely on the lattice: 
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a ~ /i^^ being the lattice cut-off, without hav- 
ing to introduce an intermediate (renormaliza- 
tion) scale parameter. This approach automati- 
cally takes care of all mixing effects, and it avoids 
renormalon ambiguities altogether. 
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Figure 1. The eigenvalues w„ against n. 



Figure 2. The residual error R^ against n. 



In this talk we shall present first results of an 
all-lattice calculation of higher-twist corrections 
to the structure function of the nuclcon. 

2. WILSON COEFFICIENTS 

Let us concentrate on the Wilson coefficients 
first. The basic idea [ || is to compute W^,y, and 
the operator matrix elements in the OPE going 
with it, for off-shell quark states in Landau gauge, 
and to extract the Wilson coefficients from this in- 
formation. Because the quarks are off-shell (and 
euclidean), we are far from the Bjorken limit, so 
that higher-twist operators are not suppressed. 

We consider a system of quark states with mo- 
menta Prm 'fn = 1, • • • , M. The photon momen- 
tum q is held fixed. We drop the Lorentz in- 
dices, and we label the operators and the as- 
sociated Wilson coefficients by n = 1, •••,]¥. 
The matrix elements of the operators between 
states of momentum pm are denoted by OJJ'", 
and the corresponding current matrix elements 
are called VF^". Note that both, 0^- and W^^^ 
are 4x4 matrices. The problem is then to solve 
the N X (Af X 16) system of equations 
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for the Wilson coefficients C„. In the following 
we take into account all quark-bilinear operators 
with up to three covariant derivatives. 



Generally, not all operators are independent, so 
that the system is overdetermined. Let us write 
(0)as 



OC = W. 



(3) 



To solve (|3|), we apply a singular value decompo- 
sition, which is the standard method for solving 
overdetermined sets of equations [ Q . Writing 

= UloV'^, 

where U {V) is a column-orthonormal (16 x M) x 
N {N X N) matrix, and oj — diag(a;„) with pos- 
itive real eigenvalues u;„ arranged in descending 
order, the solution to (^ is 

C = Vd[(ig{l/uJn)U^W. 

In a calculation with infinite precision one would 
find Un = for any n exceeding the number of 
independent operators, and in these instances one 
would replace l/cj„ by zero. In a calculation like 
ours we can at most hope for a sharp drop of the 
(jj's down to the level of the noise. 

In total there are 1360 operators to be consid- 
ered, which by symmetry reduce to 93. We have 
considered 70 different momenta, thus resulting 
in a 93 X 1120 system of equations. The calcu- 
lations are done on a 32^ lattice at (3 = 6 using 
Wilson fermions. In Fig. 1 we show a typical pic- 
ture of the eigenvalues w„ . We see a sharp drop 
at n ~ 63. In Fig. 2 we show the residual error 
R^ = \W ~ OC\^/\W\^ as a function of n with 
l/cjn+i, • • • , l/uiN set to zero. We see that the 



error does not decrease significantly anymore for 
n ^ 63, indicating that we have 'hit the noise', and 
that we may truncate the system at this value of 
n, what we will do. 



3. STRUCTURE FUNCTION RESULTS 

The next step is to compute the nucleon ma- 
trix elements of the operators entering the OPE 
(M. This, combined with the Wilson coefficients 
evaluated in the last section, then gives us the 
desired hadronic tensor VF^i/ of the nucleon. For 
this calculation it is sufficient to consider a few 
different nucleon momenta only. 



o c^/q' 



From W^jy we can derive moments of the nu- 
cleon structure functions. Here we shall restrict 
ourselves to the lowest non-trivial moment of the 
unpolarized structure function. Projection onto 
spin two states gives [ H] 

jde sin^e C] (cos6')W^AA oc jdx {2xFi-Fl/2), (4) 

where C2(cos0) is a Gegenbauer polynomial, and 
cos 9 = pq/\p\\q\. Similar results are obtained for 
other combinations of /x, v. 

The nucleon matrix elements are computed on 
a smaller, 16^32 lattice (at the quoted value of 
/3), and the combined results are extrapolated to 
the chiral limit. The calculations are done in the 
quenched approximation, so that our numbers re- 
fer to non-singlet structure functions. 

We denote the r.h.s. of (||) by A^2(9^)- In the 
language of the parton model this corresponds to 
the moment (x). We write 

M2{q^) ^m2 + C2/q^, 

separating the power correction from the leading, 
logarithmically varying contribution. This makes 
sense, because we can identify those higher-twist 
contributions which are logarithmically behaved, 
and those which arc power behaved. 

We present our results in Fig. 3. All our num- 
bers should be regarded as preliminary. The first 
observation we make is that power corrections are 
large. They are positive, and at g^ = 5GeV^ 
they amount to « 30% of the total contribution. 
The second, and equally important observation is 
that there is strong mixing between operators of 
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Figure 3. A^2('Z^) for the proton as a function of 
q^ . The hatched strip shows our old result [ ||] 
(where the width of the strip indicates the error) 
without higher-twist effects included. 



twist two and four. This effect reduces the leading 
contribution (m2) substantially, if compared with 
the result of our previous calculation [ |^ based on 
perturbative Wilson coefficients and renormaliza- 
tion constants. We now find much better agree- 
ment of this quantity with the phenomenological 
valence quark distribution functions. 
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